Large-scale prediction of the parity distribution in the nuclear level density and 

application to astrophysical reaction rates 



D. Mocelj/'B T. Rauscher/'Q G. Martmcz-Pincdo, 2 K. Langanke, 2 ' 3 
L. Pacearescu, 4 '0 A. Faessler, 4 F.-K. Thielemann, 1 and Y. Alhassid 5 

1 Departernent fur Physik und Astronomie, Universitat Basel, Basel, Switzerland 

2 GSI Darmstadt, Darmstadt, Germany 
3 Institut fur Kernphysik, Technische Universitat Darmstadt, Darmstadt, Germany 
J 'Institut fur Theoretische Physik, Universitat Tubingen, Tubingen, Germany 
5 Center for Theoretical Physics, Sloane Physics Laboratory, Yale University, New Haven, CT, USA 

A generalized method to calculate the excitation-energy dependent parity ratio in the nuclear 
level density is presented, using the assumption of Poisson distributed independent quasi particles 
combined with BCS occupation numbers. It is found that it is crucial to employ a sufficiently large 
model space to allow excitations both from low-lying shells and to higher shells beyond a single 
major shell. Parity ratios are only found to equilibrate above at least 5-10 MeV of excitation en- 
ergy. Furthermore, an overshooting effect close to major shells is found where the parity opposite to 
the ground state parity may dominate across a range of several MeV before the parity ratio finally 
equilibrates. The method is suited for large-scale calculations as needed, for example, in astrophys- 
ical applications. Parity distributions were computed for all nuclei from the proton dripline to the 
neutron dripline and from Ne up to Bi. These results were then used to recalculate astrophysical 
reaction rates in a Hauser-Feshbach statistical model. Although certain transitions can be consid- 
erably enhanced or suppressed, the impact on astrophysically relevant reactions remains limited, 
mainly due to the thermal population of target states in stellar reaction rates. 
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I. INTRODUCTION 

Knowledge of the nuclear level density in general and, 
more specifically, at low excitation energies is of interest 
for a number of reasons. Comparison to experimental 
level densities helps testing nuclear structure models [l[ . 
In reaction theory, the nuclear level density is an impor- 
tant ingredient both for the determination of the relevant 
reaction mechanism and for the calculation of reaction 
cross sections [2[ . Because of the low effective interaction 
energies in astrophysical applications, the level density at 
low excitation energy is usually assumed to be crucial to 
determine astrophysical reaction rates. 

The level density p as a function of spin J, parity 7r, 
and excitation energy E can be written as 

p(E,J,7r)=V(E,ir)T(E,J) Ptot (E) , (1) 

with the spin projection T . In most previous applications 
to astrophysics (e.g. (H), equally distributed parities with 
P(E,tt = ±1) = 1/2 for both even and odd parity have 
been assumed at all energies but it is obvious that this 
assumption is not valid at low excitation energies. This 
work focuses on the determination of the parity projec- 
tion factor V and its implication on the calculation of 
reaction rates for astrophysics. 
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In astrophysical applications usually different aspects 
are emphasized than in pure nuclear physics investiga- 
tions. Explosive burning produces a large number of un- 
stable nuclei for which experimental information is un- 
available. Thus, the study of such models requires pre- 
diction of nuclear properties for a large number of nuclei, 
several thousands when considering all nucleosynthesis 
processes. On the other hand, reaction rates are obtained 
by energy-averaging cross sections in the relevant energy 
range (see Eq. [2^)) . given by the so-called Gamow window 
p|. Therefore, small deviations from the true cross sec- 
tions may cancel out. This also may dampen the impact 
of variations in the predicted nuclear properties, such as 
the nuclear level density. 

In principle, the nuclear level density should be ex- 
tracted from microscopic models. However, conventional 
shell model calculations of level density are limited to 
the mass range A < 50 0, i, @] because of the combina- 
torial increase of the dimension of the model space with 
the number of single-particle levels and/or the number of 
valence nucleons. Large-scale shell model calculations of 
level densities are possible in the framework of the shell 
model Monte Carlo (SMMC) method @, 1, @. Most 
SMMC calculations have been carried out in one full ma- 
jor shell. However, they can be extended to higher ex- 
citation energies by including all other shells within a 
mean- field approximation [l[. The SMMC calculations 
are in general time-consuming and are difficult to carry 
out for the large number of nuclei required in large-scale 
astrophysical applications so far. Although computers 
have become faster, a consistent microscopic description 
of all required properties for all nuclei is still not feasi- 
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ble. This is, of course, also due to the still insufficient 
knowledge of nuclear interactions in general and its ef- 
fects in neutron- or proton-rich nuclei. Therefore, many 
approaches for calculating astrophysical reaction rates 
make use of more phenomenological treatments. This 
does not imply that they are just mathematical fits to a 
given property, made in the region of stability and then 
extrapolated to unstable nuclei. Instead, they are based 
on some physical insight in the origin of the given prop- 
erty which makes it possible to extend its application 
also to unknown nuclei with only few parameters to be 
adapted. Among the traditional examples of such ap- 
proaches is the shifted Fermi-gas model [l(J [H| °f the 
nuclear level density and its variations. As long as nu- 
cleons can be assumed to exhibit a Fermi-gas behavior, 
the energy dependence of the level density can be de- 
scribed by a few parameters. This has been proven re- 
peatedly by comparison to experimental data but also 
in comparison to microscopic calculations [l|, [H, [H, [l4| • 
Also the behavior of its parameters can be determined 
by applying known physical facts and/or by extracting 
them from different types of calculations 0, [l5| ■ Another 
example of the combined phenomenological approach is 
the macroscopic-microscopic Finite Range Droplet Model 
(FRDM) [r| which still proves to be the most successful 
model to predict nuclear masses. 

In the spirit of the above we present a macroscopic- 
microscopic method to determine the parity ratio of nu- 
clear levels from low to high excitation energy which can 
be easily applied to a large number of nuclei. Providing 
the factor V as a function of excitation energy - as we do 
in the following - has the additional advantage that it can 
be combined with any total level density p t ot from any 
other approach to determine the density of odd and even 
parity states. The method is based on a quasi-particle 
model proposed in Ref. [l7| which we extend by allowing 
odd particle numbers, inclusion of all shells up to the 11 
hio oscillator shell, and including excitations between all 
considered shells. Preliminary results were reported in 

[H d, us mi . 

The method is introduced in Sec. [Til with a discus- 
sion of some interesting effects found in the parity ratios 
included in Sec. Ill Dl Sec. IIIII then focusses on the as- 
trophysical relevance of the parity distribution. Astro- 
physical rates computed with the new parity ratio are 
compared to standard rates widely used in astrophysical 
applications in Sec. IIII CI The paper is concluded with a 
summary in Sec. IIVI 



sumed to be binomial [l3[. The probability to occupy 
n out of r levels is therefore given by: 



B(r;n)= ( )p n (l - p)" 



(2) 



where the probability of occupying a level is denoted by p. 
We follow [T?! by replacing the binomial distribution by 
a Poisson distribution P(s). This approximation holds 
provided the number of levels r is large and the prob- 
ability p is small, with the product rp finite. This can 
always be achieved by clustering the single particles into 
two groups according to their parity and counting only 
particles in the group having the opposite parity to the 
last occupied level. The probability to find n particles in 
that group is then 



Pin) = 



-f 



(3) 



where the average number of particles in the group is 
given by /. Obviously, neutrons and protons have to be 
treated separately. For an even number of nucleons the 
probability to find the whole system in a positive parity 
state is therefore given by 

P+= J2 J —e- f = cos\xfe-f, (4) 



and to find the system in a negative parity state by 

P- = V i- e -f = smhfe- f . (5) 

n,odd 

The probabilities can be related to the total partition 
function 



z = z + + z~ 



(6) 



and the partition functions for odd and even parity 
states, Z~ and Z + , by 



P+ = 
P- = 



Z+ 

~z~' 

z~ 

~z' 



(7) 



This leads to the expression used in [17| for even-even 
nuclei: 



II. DETERMINATION OF PARITY RATIOS 



Z+ 



tanh /. 



(8) 



A. Basic approach 

We start from the assumption of statistically indepen- 
dent particles at finite temperature. Because single par- 
ticle (s.p.) levels can either be occupied or empty, the 
probability distribution for the occupation can be as- 



It can be shown [22| that it can be completely general- 
ized for even and odd numbers of particles by denoting 
the total ground state parity of a nucleus by g and the 
opposite parity by s: 



ps ^ Z s 

~pl ~ ~Z~9 



tanh /', 



(9) 



3 



where /' = f a + f p is now computed from the sum of the 
individual average particle numbers / of neutrons and 
protons. 

In order to determine the desired parity factor V the 
ratio of the level densities p g , p s with different parity has 
to be known. Applying the well-known Laplace trans- 
form of the partition function and employing the saddle 
point approximation, the ratio for an excitation energy 
E is given by [22| 



f)s_ _ (3 s Z s C 9 j0 S _ l3 a )E 

Pg 139 Z9 V C° 
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with the heat capacities C 9 , C s . In general, because 
of the different number of particles in each group, the 
inverse nuclear temperatures /3 s , (3 9 will be different for 
the same excitation energy E. The required partition 
functions are determined by 



Z 9 = 

z s = 



1 



1 + tanh/'' 
1 
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(11) 



tanh /' 



The thermal energy E and the heat capacities can be 
derived from the standard thermodynamic relations 



E = - 
C = -f3 2 



d\uZ 



dft 

d 2 InZ 

dp 2 



(12) 



Finally, the parity projection factor V is given by 

1 



V{E 1 TT = K g )=V g = 

V{E,ir = ir s )=T s = 



Pg 

Ptot 
Ps 
Ptot 



i + e 
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with 



Pg 



(13) 
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and -Kg being the ground state parity of the nucleus while 
ir s being the opposite parity. 



B. Determination of the mean occupation number 

The central quantities / n , f p were determined sepa- 
rately for neutrons and protons. In the following, we 
denote them by / for simplicity and imply that the 
two kinds of nucleons were treated separately with their 
respective s.p. levels and particle numbers. Assuming 
Fermi-Dirac distributed particles, the occupancy f k of 
each s.p. level with energy e k given by 
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(15) 



and the average number of particles by 



= £/* 



(16) 



Proceeding as before by dividing the s.p. levels into a 
group exhibiting the same parity ir g as the last occupied 
state below the chemical potential p, and another group 
with opposite parity ir s , the mean value / in the Poisson 
distribution is then given by 
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(17) 



The chemical potential is found iteratively from the par- 
ticle number equation 
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With decreasing nuclear temperature T nuc pairing in- 
teractions become increasingly important in nuclear sys- 
tems, leading to nucleon pairs in the ground state and at 
low excitation energies. The breaking of such pairs re- 
quires additional energy. Hence, the Fermi-Dirac distri- 
bution will not be able to describe the occupation prop- 
erly at high values of j3 = 1/T nuc . For example, this 
can be seen in Fig. 3 of [ItJ where a comparison be- 
tween Fermi-Dirac and BCS occupancies and resulting 
partition ratios Z_/Z + are shown for even-even nuclei. 
At high nuclear temperatures the microscopic distribu- 
tions are well described by Fermi-Dirac statistics but for 
lower temperatures deviations start to appear. For those 
cases of even-even nuclei, the probability to find an even 
number of particles in the n s parity group is clearly en- 
hanced, odd numbers of particles are suppressed as seen 
in the comparison to the SMMC result in [T3|. Similar 
effects were reported in [l[. 

Pairing effects can be included by using the well-known 
Bardeen-Cooper-Schrieffer (BCS) formalism (23j, intro- 
ducing the BCS occupancy for quasi-particles 



/f cs = 1/(1 + eMPE k )) 



with the quasi-particle energy 



E k = v/(e fe - pf + A 2 



(19) 



(20) 



The chemical potential jl = p(P) and the pairing gap 
A = A(/3) are determined for each nuclear temperature 
?nuc = 1//3 by solving a non-linear system of equations, 
the well-known particle number and gap equation 



2^V E * \ 2 



2 
A 



2^ — tanh ( ^E k 



2^E k 



(21) 



where G is the usual effective pairing coupling constant 
which determines the zero-temperature pairing gap A. 
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Although the single particles are no longer independent 
in the paired solution, the Poisson approach of Sec. Ill Al 
can be kept, assuming that the quasi particles are sta- 
tistically independent here. For the calculation of the 
mean occupation in the tt s parity group, it is sufficient to 
consider only contributions from quasi-particles, as con- 
densed pairs do not contribute to any parity change of 
the system. This leads to 



/ 



fcG7r 3 



fc67T 3 



exp(l3E k ) 



(22) 



It should be noted that here the chemical potential p, and 
the pairing gap A are temperature-dependent, i.e. have to 
be known as a function of f3 or T nuc . The BCS equations 
are solved iteratively for each inverse temperature f3. 

At sufficiently high nuclear temperature, pairing be- 
comes negligible and there is a phase transition from the 
BCS to the Fermi-Dirac regime. However, this transi- 
tion is not smooth when only using the above description 
because the BCS equations break down above a critical 
temperature T c — l/(3 c . Due to the phase transition 
occuring at T c , the derivatives of the partition function 
Z (containing the mean occupancy /) are discontinuous 
when moving from the BCS regime to the Fermi-Dirac 
regime. This induces numerical problems in the solution 
of Eq. ([TO)) via Eq. (| 12[) . especially also because the dif- 
ferent parity groups exhibit different temperatures. To 
avoid the problematic temperature range would be too 
limiting because the range of excitation energies most 
interesting for studying the parity ratios would be ex- 
cluded. Therefore, we employed an approximation by 
extrapolating the pairing gap to higher nuclear tempera- 
tures. In fact, the transition from the BCS to the Fermi- 
Dirac regime should be smooth in nuclei because the pair 
coherence length is much larger than the nuclear radius 
In such systems with dimensions smaller than the 
coherence length, fluctuations in the order parameters 
become important and wash out the discontinuity. 

As the pairing effects vanish for high nuclear temper- 
atures, the pairing gap itself should approach a value of 
Zero, thus recovering the Fermi-Dirac shape of the distri- 
bution (see, e.g., Ref. (2f| for experimental indications of 
this effect). This behavior of the gap could be rigorously 
obtained, for instance, by using the static path approx- 
imation, integrating over the static fluctuations of the 
gap [26|. However, as has been shown in SMMC calcula- 
tions [3, [2^ | , this behavior can be simply parameterized 
in terms of a Fermi-type dependence: 



A ht {P) = A ht (T, 



An 



l + cxp(-(T n 



(23) 

where An. denotes the pairing gap for T nuc = 0. The pa- 
rameters T^ uc and a were fitted so that the function ap- 
proximated as closely as possible the behavior of the BCS 
gap below the critical temperature. As starting values we 
chose T/ mc = 0.8T C and a = (2^-T nnc )/ln(A /A£F[- 
1) which yielded a good fit to the BCS pairing gap 
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FIG. 1: Dependence on nuclear temperature of the proton 
pairing gap A in 66 Zn and its fit by a Fermi- type function 
(see text). 



A BCS = A below jv^ Thig resultg 

in a suppression of 
the pairing gap across a temperature range of the order 
of MeV, as also previously found in 0, E2H3] • For illus- 
tration, Fig. [U shows a comparison of the BCS solution 
and the fit for 66 Zn. The results are not sensitive to small 
variations in the choice of the starting values. These fits 
were subsequently used in the calculation of the chemical 
potential fi and finally of the mean occupation number / 
by applying Eq. (|2"2"|) for all values of (3. 



C. Input quantities and consistency 

The parity projected partition functions Z 9 , Z s are re- 
lated to each other, to the total partition function Z, and 
to the average occupancy / by Eqs. © and ©. Thus, 
the total partition function Z and the average occupation 
number / remain to be determined. 

The determination of / as described in the previ- 
ous section implicitly requires knowledge of the pairing 
strength G and the s.p. levels. For systems with an even 
number of particles, the starting value Ao of the pair- 
ing gap was extracted from odd-even mass differences as 
described in Q. This value was then used to determine 
the effective pairing strength G by solving the coupled 
Eqs. d2JJ for T nuc = 0. Further application of Eqs. (j2Tj) 
with fixed G at finite nuclear temperatures yields the 
temperature-dependence of A and the required occupa- 
tion of quasi-particle states. Differently from the treat- 
ment of the pairing gap Ao in Q , here we set Ao = for 
systems with odd numbers of particles because the pair- 
ing effects are expected to be weak and quickly vanishing 
with increasing T nuc . This means that a pure Fermi- 
Dirac distribution is used for odd particle numbers with- 
out the need to fit the phase transition. 

The s.p. levels were calculated in a deformed Saxon- 
Woods potential [28| with parameters from [2^| which 
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FIG. 2: Ratio of the level densities with odd and even parities 
as a function of excitation energy _E X in MeV in the nucleus 
56 Fe, calculated with different sizes of the model space. It is 
evident that it is not sufficient to only include the pf + gg/2~ 
shell. 



reproduce experimental data well [30L l3l| . The same de- 
formation was used as in 0, [32l . |33|, taken from . 

Fig. [2] illustrates the dependence of the parity distribu- 
tion on the single-particle model space in 56 Fe. For this 
nucleus, the addition of the 2d 5 / 2 orbital to the pf + gg/2~ 
shell leads to a deviation over 10% in the parity ratio 
above excitation energies E x of about 15 MeV. Further 
inclusion of higher states such as 1 177/2 does not change 
the result for 56 Fe (the lines for the results overlap in the 
figure) but will be important for heavier nuclei. There- 
fore we had to include all major shells up to 11 fwj. 

We also see in Fig. [5] that levels usually considered 
as inert core can also have an important effect. This is 
clearly seen in the result accounting for excitations across 
the full sd + pf + gds-shells. Here we find deviations of 
at least 10% already at E x ~ 7 MeV and larger devia- 
tions at higher excitation energies. This is mainly due to 
excitations from the sd-shell. 

As a further ingredient, a total nuclear level density 
Ptot is required in the determination of the total partition 
function Z — Z((3, ptot)- In principle, any total level 
density determined in any approach could be supplied 
here. Since we want to apply the derived parity factor 
to the level density of [H, we took the total level density 
and the level density parameters from there. 

In a fully consistent model, s.p. levels determine the 
level density assuming all effects have been properly ac- 
counted for and provided that the level density can be 
computed, e.g. in a shell model, for all nuclei at all 
relevant excitation energies. In principle, this can be 
achieved, e.g., within the approach of Ref. [l|, but is dif- 
ficult to carry out for a large number of nuclei. For prac- 
tical purposes it is therefore necessary to include level 
densities derived in other approaches, such as the shifted 
Fermi-gas. Proceeding in this manner, however, great 
care has to be exercised to achieve consistency between 
the s.p. structure, in our case obtained from the Saxon- 
Woods model, and the input level density. Inconsistency 
will lead to spurious effects in the parity ratio as shown 
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FIG. 3: Influence of a variation of the input level density 
Pinput on the calculated parity ratios in 70 Zn. The full line 
is obtained by using the standard level density p Bt d 0. The 
dashed and dotted lines correspond to variations of the input 
level density by factors 2 and 5, respectively. The s.p. levels 
are held fixed. 



in Fig. GS 

Realistically one has to expect that both the s.p. level 
structure from the Saxon- Woods potential and our to- 
tal level density bear inherent uncertainties. Thus, when 
striving for consistency it is anticipated that both inputs 
should be adapted. Here, we used an iteration method 
which only varies the level density and not the Saxon- 
Woods potential, assuming that the s.p. structure is cor- 
rect. Although this does not limit the applicability of our 
results, it has to be emphasized that the variation of the 
level density pinput = Ptot shown in Fig. [3] and found in 
our iteration procedure is not due to the uncertainty of 
Ptot alone but rather is supposed to contain the combined 
uncertainties in p to t and the s.p. levels. 

In each iteration step the parity projected level- 
densities 

p g (E) = -l=^eM(3 9 E + \nZ^ , 

V IttC 3 

Ps (E) = -l=p s exp(p s E + lnZ s ) (24) 

V 2ttC s 

were calculated with an input level density p; n put- The 
total level density in this step is then given by the sum 
of the parity projections: 

p^ =p g + p s . (25) 

In the first iteration step the standard level density from 
Ref. [2| was used for pi„ pu t- In each following step, p^ of 
the previous step became the new pi npu t and the proce- 
dure continued until p^ and p; np ut converged. Conver- 
gence was usually achieved quickly, typically in the third 
iteration step. 
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FIG. 4: Ratio of the parity projected level densities as a func- 
tion of excitation energy for Ni isotopes around the neutron 
sd-shell closure. 
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FIG. 5: Ratio of the parity projected level densities as a func- 
tion of excitation energy for Ni isotopes around the neutron 
p/-shell closure. 



D. Results and discussion of selected examples 

Using the approach described in the previous sections, 
we have calculated ratios of parity-projected level den- 
sities for all nuclei from Ne to Bi and from proton 
dripline to neutron dripline, adoptin g th e large HTia; 
model space. All target nuclei in [32[ l33j] are covered 1 . 
The results can be obtained from the American Institute 
of Physics' EPAPS [34j or directly from the authors at 
|http : / / nucastro . org) 

Representative for our results are Figs. I3HH] It can 
immediately be noticed that the parities equilibrate only 
at excitation energies above about 5—10 MeV or even 



1 It should be noted that this also impacts, for example, the re- 
sulting driplines as the same inputs as in these references were 
used here. 



higher (see, e.g., Fig. [4|, even for the more heavy nuclei. 
This underlines the importance of using a parity projec- 
tion factor V 7^ 1/2 at most energies of astrophysical 
interest. 

Surprising at first glance is the effect we found at closed 
shells. The parity ratio can overshoot the equilibrated 
value of Ps/Pg = 1 by 20-60 % before finally reaching 
the equilibrium value at higher excitation energies. The 
overshooting can be understood by the fact that parity 
can only be changed by excitations between s.p. levels of 
different parity. When an oscillator shell is completely 
filled, i.e. the next available subshell exhibits opposite 
parity, any excitation will change the parity of the sys- 
tem, resulting in a dominance of opposite parity states. 

While studying the following examples, it should be 
kept in mind that both neutron and proton excitations 
contribute to the parity ratio. When comparing the ra- 
tios within an isotopic chain, however, the proton con- 
tribution will remain the same and differences between 
isotopes can be attributed to the change in neutron num- 
ber. Similar considerations apply to the comparison of 
isotones and the changing proton number. For example, 
Fig. 2] shows the evolution of the parity ratio in Ni iso- 
topes around the neutron sd-shell closure. The sd-shell 
is filled completely with neutrons for 48 Ni. Every excita- 
tion of 48 Ni populates a single particle level of opposite 
parity leading to a maximal parity change. The forma- 
tion of the peak can already be seen for 47 Ni but the 
position of the peak is shifted towards smaller excitation 
energies. This energy shift is about 1 MeV reflecting the 
energy that is needed to break a neutron pair prior to 
excitation. The isotopes 49 Ni and 50 Ni, which both al- 
ready populate the p/-shcll in neutrons, equilibrate at 
much higher excitation energies. Parity change can be 
achieved by excitations to the g$ / 2 shell or by excitations 
from the sd to the pf shell, requiring higher energy on 
average than for the preceding isotopes. 

The evolution of the parity ratio for Ni isotopes at 
the interface of the pf and gds shell is shown in Fig. 
[5] Again, a similar behavior as described above can be 
observed. The N = 40 neutron shell is completely filled 
for 68 Ni. Each excitation from the last occupied 2p 1 / 2 
level with negative parity will populate levels from the 
gds shell with positive parity resulting in a parity-change 
of maximal amplitude. The formation of the peak can 
already be seen for 66 Ni. The position of the peak is 
shifted to higher energies as a larger gap between the 
2/5/2 level and the 59/2-level has to be bridged. 

The transitions between the same shells as above in 
Fe isotopes is shown in Fig. [5] Moving from 56 Fe, where 
the neutron pj-shell is only half-filled, to 66 Fe, where the 
2p 1 / 2 -shell is completely filled, it can be seen that the 
ratio approaches unity for lower values of the excitation 
energy as one approaches the N = 40 shell closure. As 
the parity can only be changed by excitations either from 
the sd to the pdf shell or from the pf to the gds shell, 
the ratio will equilibrate faster with increasing neutron 
number as the gap between the last occupied orbit in 
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FIG. 6: Ratio of the parity projected level densities as a func- 
tion of excitation energy for Fe isotopes up to N = 40. 
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FIG. 7: Evolution of the ratios of the parity projected level 
densities as a function of excitation energy for Sr isotopes 
around the neutron p/-shell closure. 
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FIG. 8: Evolution of the ratios of the parity projected level 
densities as a function of excitation energy for Sn isotopes 
around the neutron ods-shell closure. 



5 10 15 20 5 10 15 20 



- ?2 Mo 


I 




I 



5 10 15 20 5 10 15 20 

E x (MeV) 

FIG. 9: Evolution of the odd- to even-parity ratio within the 
N = 50 isotones in the vicinity of the proton p/-shell closure. 



the p/-shell and the gds-shell will decrease. For 66 Fe, 
where the pf shell is completely filled, a pronounced peak 
around 8 MeV occurs. 

A similar behaviour is illustrated by Fig. [7] for the evo- 
lution of the parity ratio in Sr isotopes at the pf and gds 
interface. Again, the ratio reaches a maximum in 78 Sr 
with the filling of the neutron p/-shell. A comparison to 
the Fe isotopes with the same neutron numbers in Fig. 
P shows the impact of deformation. The FRDM model 
[10 ] predicts the Sr isotopes shown here to be strongly 
deformed. The resulting level splitting leads to a wider, 
less pronounced peak in the Sr isotopes as compared to 
the Fe isotopes. 

Another interesting case is seen in Fig. [5] which dis- 
plays the evolution of the parity ratio in Sn isotopes in 
the vicinity of the gds- and pfh-shell. The /i n / 2 -subshell 
is half filled at 120 Sn due to the fact that it is located 
below the 2d 3 / 2 and Ss^-subshells, the latter two be- 
longing to the gds-sheW. Comparatively small excitation 



energy is required to change the parity by either moving 
neutrons from the /iii/2-subshell to the higher ds-states 
or by populating the /ii^^-subshell with neutrons from 
lower gd-states. The effect of the latter is largest when 
the /111/2-subshell is least occupied, the effect of the for- 
mer becomes larger with increased occupation. That is 
why there is an overshooting peak found in all Sn iso- 
topes shown in Fig. [5J with the one in 120 Sn being the 
most pronounced because both excitations contribute al- 
most equally. However, due to the mixing of gds- and 
/i n /2-states, the overshooting is not as strong as the one 
found for 66 Fe at the ^/-closure. 

Above, the evolution of the parity ratio was discussed 
in isotopic chains, where each parity change results from 
an excitation of neutrons only. A similar effect is found 
for protons. For instance, Fig. [5] shows the evolution of 
the N = 50 isotone parity ratio. The proton numbers for 
the four shown nuclei (Y, Zr, Nb, and Mo) are 39, 40, 
41, and 42, respectively. Again, a peak can be observed 
in 90 Zr as the proton p/-shell is completely filled for this 



8 



nucleus. Weak deformation is predicted for these nuclei 
which leads to a suppression of the peak. 



III. APPLICATION IN ASTROPHYSICS 

A. Introduction 

In the following we discuss the implications of a re- 
alistic parity distribution in astrophysical applications. 
The astrophysical reaction rate is the central quantity in 
the nuclear reaction networks employed to follow nucle- 
osynthesis in different astrophysical environments. For 
nucleon-nucleus and nucleus-nucleus reactions it is given 
by folding the reaction cross section with the Maxwell- 
Boltzmann velocity distribution of the interacting nuclei 



nM(kT*y 



- OO 



(26) 

The stellar plasma temperature is denoted by T*, the 
reduced mass of the interacting nuclides by M, and their 
number densities by n\, n 2 . The stellar cross section is 
defined as 



<7*(£ c ,n.) 



E M (2J^ + 1) exp(-E^/kT*) £„ a^(E c . m .) 
E„(2J^ + l)exp(-^/fcT* 



(27) 

including sums over excited states /i, v in the target and 
final nucleus, respectively, with J M and denoting spin 
and excitation energy of the //-th target state. 



B. Relevant reaction theory 

For the majority of reactions relevant in astrophysics, 
the reaction cross sections a^ v have to be calculated with 
different reaction models as experimental information of- 
ten is not available in the relevant energy range or even 
impossible to obtain due to the short half-lives of the 
involved nuclei. Even if data is available, it can only 
provide a measure of cr^ = ^ v tr A " y and thus does not 
provide the required a*. Depending on the number of 
resonances in the Gamow window [2|, one has to con- 
sider direct reactions, the influence of few resonances, or 
an average over many resonances (35J. Here, we want to 
focus on the latter and on direct capture. 

Provided there is a high number of resonances at 
the relevant interaction energy, the formation and de- 
cay of a compound nucleus can be described by aver- 
aged transmission coefficients in the statistical model or 
Hauser-Feshbach formalism (36j. The transmission coef- 
ficients describing the decay of a compound state k via 
a process y - with y usually being emission of nucle- 
ons, a particles, or 7 quanta - include sums of the type 
k T j k -n k _ j2 v k T J k * h ^J"ir" _ Depending on the nucleus, 




A A 



FIG. 10: Energy scheme of a compound nucleus reaction 
A + pro — » C* — > ej + B. For a capture reaction on tar- 
get nucleus A the ejectile "ej" will be 7 radiation and the 
final nucleus B will be identical to the compound nucleus 
C. The standard reaction Q value is the difference be- 
tween the separation energies of the projectile and the ejectile 
Q = Spro — S c j. When considering individual transitions be- 
tween excited states in the target and the final nucleus the 
released energy becomes Q = (S pTO + — (S c j + E B ), with 
the excitation energies E^, E B . The excitation energy of the 



accessed compound state is then E c 



Spro ~\~ Ep ro ~\~ Ea 



with E pro being the projectile energy. The decay of the com- 
pound state k is computed by a sum of transitions to indi- 
vidual states and by integrating over an average level density 
' (shaded area) above the last included state. The target A can 
become excited by thermal population of the excited states 
in an astrophysical plasma. 



reaction Q value, and interaction energy -E c .m., the prop- 
erties (spin J u , parity 7r" ', excitation energy E v ) of the 
accessible states may not be known above a certain en- 
ergy E c . Then the sum is replaced by an integration over 
a level density 



I J- Tf^ J "«"(Q + E c 



xp{E v ,J v ,-K v )dE 



E v ) x 
(28) 



The integration starts from E c = in most cases far 
from stability where no excited states are known. The 
level density p(E h ' , J v , n v ) includes the assumption of the 
parity distribution. It should be noted that this is the 
level density in the final nucleus which is identical to the 
compound nucleus only for capture reactions. The inter- 
action energies -E c .m. in astrophysics typically are smaller 
than 1 MeV for neutrons and 5-10 MeV for charged pro- 
jectiles. A comparison with our results from the previous 
section shows that the parity ratio obviously has not equi- 
librated at the relevant excitation energies, even though 
the reaction Q value has to be added to these energies. 

Considering thermal excitation of the target states in 
the transmission coefficients for the formation of the com- 
pound nucleus in state k a relation similar to Eq. 
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be derived for the transmission coefficients T y i connect- 
ing target and compound state via a process y' , this time 
inserting the level density in the target nucleus. However, 
due to the exponential suppression of the population of 
the excited target states (see Eq. [27f significant contri- 
butions to the integral will only stem from the lowest 
energies. 

The central assumption of the Hauser-Feshbach model 
is that the compound state at energy E k is highly de- 
generate. Due to the large number of overlapping reso- 
nances, the excited compound nucleus can have any com- 
bination of spin and parity. The possible combinations 
will be weighted by the transmission coefficients k T y i de- 
scribing the formation process and, e.g., accounting for 
spin selection rules. However, in principle one has to 
sum over all compound spins J k and parities n k when 
calculating the statistical model cross section 

krp,J k TT k krp,] k lT k 

oc £ (2J* + 1) k *j^ . (29) 

z=y,y 

The denominator contains a total transmission coefficient 
for each state J k ir k , including all possible deexcitation 
processes, e.g. emission of protons, neutrons, a particles, 
and 7 quanta. The scheme of a compound nucleus reac- 
tion is sketched in Fig. rXUl 

Direct reactions proceed in a different manner, without 
formation of an excited compound nucleus. A nucleon or 
nucleon group is directly transferred from the projectile 
to the final state without excitation of any other nucleons 
in the system 35, 37, 38]. This mechanism will dominate 
at very high interaction energies where the compound 
nucleus formation is suppressed due to the short interac- 
tion timescale and at low energies in the absence of reso- 
nances, again suppressing the formation of a compound 
nucleus. In astrophysics, direct capture is expected to 
be important close to the driplines because of the low 
particle separation energies. The cross section a^, for 
direct capture of a projectile on a target in state fj, into a 
final state v can be calculated, e.g. in the potential model 
[35I I371 . [38| or in the Lane-Lynn formalism [3^]. Again, 
the total capture cross section is the sum of the contri- 
butions of the different transitions to discrete states 

xa^ c (Q + E c . m . -E V ,J^,^,E V ,J V ,7T V ), (30) 

with C and S denoting the isospin Clebsch-Gordan co- 
efficient and the spectroscopic factor of the final state, 
respectively. Because of the importance of direct capture 
close to the driplines only, the discrete states are usually 
experimentally unknown and have to be extracted from 
microscopic calculations. However, it has been shown 
[40| that the resulting cross section is very sensitive to 
the predicted properties of the states and that thus dif- 
ferent microscopic models yield vastly different results. 



Pa 




FIG. 11: Energy scheme of a direct capture reaction ^4+pro — > 
B + 7. The energy of the emitted 7 quantum is the energy 
difference between initial target state and final state. For 
instance, Qoo is the released energy for the ground state to 
ground state transition and is equal to the binding energy of 
the projectile in the final nucleus B. As in the case of the com- 
pound mechanism, it is possible to include an integration over 
a level density (shaded region) when the exact properties of 
the nuclear states are not known. For simplicity, only transi- 
tions originating from the ground state of target A are drawn 
here. For stellar reaction rates it would be important to also 
include transitions from thermally excited target states. 

The sensitivity is much higher than in the case of com- 
pound reactions because no averaged quantities are used. 
In order to circumvent the problem of the exact predic- 
tion of the final states it has been suggested [U [43j 
to employ averaged properties also for direct captures, 
i.e. to replace the sum over discrete final states by an 
integration over a level density 

Q 

J "- 7r " 

x a^^ JV (Q + E c . m _-E^dE v (31) 

and to use averaged quantities C and S which are 
assumed as being constant [43[ or energy-dependent 
[ill l42l . [44|. Again, the expression can be extended to 
include sums over thermally populated target states /x. 
The situation in the case of a direct capture reaction is 
shown in Fig. [TT] 

The sum over spin and parity in Eq. (|3ip usually con- 
tains only few contributing terms because of the domi- 
nance of El transitions. Because of the small Q values 
the integration does not include such a large energy range 
as in the Hauser-Feshbach calculations. These two facts 
are also the reason why the impact of a non-equilibrated 
parity ratio will be much larger for direct cross sections 
than for Hauser-Feshbach ones. The Hauser-Feshbach 
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model assumes implicitly that the compound nucleus can 
be formed with any spin and parity (although not with 
equal probability for each configuration) and that the de- 
cay of the compound state by 7- or particle-emission can 
proceed in transitions to a number of final states. This 
offers a larger number of possible spin/parity combina- 
tions through which the reaction can proceed than in the 
case of a direct process. 

Summarizing, the nuclear level density is important 
in two aspects. Firstly, to determine the dominant re- 
action mechanism. Secondly, to calculate cross sections 
accounting for transitions to a number of final states. 
It should be noted that these two aspects require the 
knowledge of the level density at different excitation en- 
ergies and in different nuclei, the compound nucleus in 
the former case and the final nucleus in the latter. In the 
following section IIII CI we study the impact of a changed 
parity distribution on reaction rates calculated in the two 
reaction mechanisms described above. 



TABLE I: Enhancement and suppression of transmission co- 
efficients relative to an equal parity distribution. 





Allowed Transition 


Factor 


TV = 7T 


1 = 






Ml 




-L ^ J 
TV 7= 7T 


t = 1 






El 





arithmetic mean between the contributions of the two 
allowed transitions when taking equidistributed parities. 
With the restriction of having only one parity available, 
one of these transitions will be suppressed completely. 
Therefore, the resulting transmission coefficient will be 
smaller or larger according to which of the transitions 
became unavailable. The enhancement factor can be ex- 
pressed as 



C. Comparison to standard rates and possible 
implications 

1. Statistical model rates 

Before employing the realistic parity distribution de- 
rived in Sec. [II] in calculations of astrophysical reaction 
rates, we want to arrive at an estimate of the maxi- 
mum impact to be expected. The level density enters 
in the calculation of the transmission coefficient as given 
by Eq. (|28p. For the standard assumption of equally 
distributed parities p{E u , J v , tt u ) = p{E v ,J v )/2 for all 
parities tt v . As an extreme case we suppose that the 
parity ratio is Zero at all energies, i.e. only states with 
the same parity tt° — n 9 as the ground state occur and 
p(E", J u ,tt v ) = for tt u ^ 7T g . This also implies that 
p(E v ,J v ,TT 9 ) = p(E v ,J v ). We further assume that the 
formation coefficient fe Xy ^ dominates the total trans- 
mission coefficient in the denominator of Eq. (|29[) and 
thus the process y and its transmission coefficients de- 
termine the cross section. This can be further simplified 
by comparing laboratory rates only, i.e. rates with the 
reaction target being in the ground state. 

Two cases can be treated separately, capture reactions 
and reactions with particle emission, i.e. y = 7 and 
y = n, p,a, respectively. In the latter case, ir k de- 
termines the allowed partial waves in the transition to 
the state in the final nucleus with tt" ', even partial waves 
for 7r fe = tt 1 ' and odd partial waves otherwise. For the 
capture reaction, with ~k v being in the compound nucleus, 
the parities select the allowed electromagnetic transition, 
i.e. El when ix k 7^ tx v and Ml otherwise. For simplic- 
ity, electromagnetic transitions with higher order are ne- 
glected. For our estimate we also assume that only s- 
and p- waves contribute to the cross section significantly. 

With the above assumptions, the transmission coeffi- 
cient calculated with Eq. (|2"5|) will be comprised of an 



1 + X 

the suppression factor is found to be 

2x 



1 + X 



(32) 



(33) 



with x < 1 being the ratio of p- and s-wave transmis- 
sion coefficient or of Ml and El transmission coefficient, 
respectively, 



jpz=3 for particles, 



for radiation. 



(34) 



It can be seen immediately that the enhancement cannot 
exceed a factor of two even when x is very small. On the 
other hand, the suppression can become strong depend- 
ing on the ratio of the transmission coefficients. Since 
usually x is smaller for electromagnetic transitions than 
for different partial waves of the particles, the largest ef- 
fects are expected for capture reactions. However, the 
effect depends on the involved parities. Table Q] shows 
the possible combinations. 

An examination of all ratios across the nuclear chart 
shows that most of them stay in the range given by the 
above defined J 7 ' and T^- . However, there are several ex- 
ceptions exhibiting an enhancement of up to a factor of 5. 
This can be understood by examining the partial widths 
of the relevant transitions which are directly proportional 
to the transmission coefficients. In the derivation above 
it was assumed that the total width is dominated by one 
of the transitions also appearing in the numerator of Eq. 
(|29|) and thus cancelling with the denominator. How- 
ever, if the total transmission coefficient is dominated 
by other transitions, e.g. proton- and a-emission in the 
case of a (11,7) reaction, then the changes in both T y 
and T y i have to be multiplied, leading to an upper limit 
of a factor of 4. Considering that also the total width 
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FIG. 12: Ratio of reaction rates as a function of stellar tem- 
perature T* (in 10 9 K) for the reaction 67 Se(n,7) 68 Se (note 
the logarithmic scale); the line denoted "lab: single/equi" 
gives the ratio of laboratory rates with the rate ro calculated 
assuming that only the ground state parity can appear in a 
nucleus whereas an equal distribution of parities was assumed 
in rate r\ and neglecting any experimentally known excited 
states. Similarly, the ratio of stellar rates is plotted as "sin- 
gle/equi". The full line gives the ratio of the "new" stellar 
rate ro calculated with the parity-dependence derived in this 
work and making use of the level information as given in Table 
III of Ref. | 3 3fl and the "standard" stellar rate n as published 
in Refs. [32l. |33|. For comparison, the ratio obtained without 
consideration of experimental level information is also plotted 
as "no exp. lev." (see text). 



may be slightly changed, the factor of 5 found in some 
cases is easily accommodated. Fig. [J2] shows the case 
of 67 Se(n,7) where the neutron width of the compound 
state k is larger than the 7 width but much smaller than 
both the proton- and a widths. The situation is slightly 
different for 68 Se(n,7), shown in Fig. [T31 where the 7 
width is much larger than the neutron width but con- 
siderably smaller than the proton- and a widths. If the 
ratio would be fully dominated by one process, then 
and J 7 * should increase and decrease, respectively, with 
increasing plasma temperature because at higher excita- 
tion energy of the compound nucleus more intermediate 
levels can be reached. Thus, either more levels can be 
populated by the preferred transition, giving rise to , 
or more transitions are missing for J-K In fact, the total 
effect stems from an interplay between the stellar tem- 
perature dependences of the formation and decay widths 
(and thus of the enhancement and suppression factors) 
and their ratios. This can lead to a complicated tempera- 
ture dependence of the rate ratios with initially declining 
and finally increasing ratios or vice versa. Examples are 
shown in Figs. [Til and[T51 

The above examples are extreme cases, including a 
comparison between rates with a parity ratio of zero and 
rates with equidistributed parities. To show more pro- 
nounced effects, no experimentally known levels in the 
participating nuclei were considered and the targets were 



FIG. 13: Same as Fig.frjbut for the reaction 68 Se(n,7) 69 Se 
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FIG. 14: Same as Fig. [131 but for the reaction 68 Se(p,7) 69 Se; 
note the different scale in the figure. 
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FIG. 15: Same as Fig. [H but for the reaction 67 Se(p,7) 68 Se 
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assumed to be in the ground state. For more realistic 
astrophysical reaction rates, the targets have to be ther- 
mally excited according to the plasma temperature which 
increases the number of possible channels for the forma- 
tion of the compound nucleus. Moreover, all statistical 
model calculations consider experimentally determined 
low-lying states, employing a theoretical level density 
only above the last included level. This will further re- 
duce the effect of a modified parity distribution but only 
for nuclei close to stability for which level information 
is available. Moreover, higher partial waves and multi- 
polarities than considered in the previous estimate will 
also contribute, thereby further washing out the effects 
seen above. Figs. IT2HT51 also show the actual ratio ob- 
tained when comparing a full calculation of stellar rates 
employing the parity distribution derived in Sec. [TT] and 
including experimental information from Table III of 1 33] 
with the standard FRDM rate set published in [32|, l33j |. 



2. Direct capture 

In principle, Eqs. (j32j) and (|33|) also apply for the en- 
hancement and suppression factors, respectively, in the 
direct capture mechanism. As already mentioned in Sec. 
IIIIB1 however, the transitions allowed by the spin and 
parity selection rules are more limited than in the com- 
pound model and thus the application is more straight- 
forward. Furthermore, there is no energy-dependent 
renormalization by a total width. Finally, due to the na- 
ture of the electromagnetic multipole operators appear- 
ing in the expression for the direct capture cross section, 
El radiation is dominating by far and all other multi- 
poles are strongly suppressed (see, e.g., )• These fac- 
tors will lead to a higher sensitivity of direct captures to 
the parity ratio. 

We defer a large-scale calculation of direct capture 
rates for low capture Q values to another paper but rather 
want to illustrate the above by a simplified example here. 
Assuming neutron capture on a tta — +1 target state the 
allowed neutron partial wave £ a in the entrance channel 
is determined by the parity selection rule 

(-l)^ = (-l)Sr B . (35) 

Since El (£~ i — 1) emission is by far dominating, this can 
be further simplified to 

(-1)^=-7T B , (36) 

with itb being the parity of the final state. Therefore 
only even partial waves will be allowed for ttb = — 1 
and only odd partial waves for ttb = +1- The latter 
relation would just reverse if the target state had odd 
parity. Because of the low projectile energies relevant 
in nuclear astrophysics, only s- and p- waves will have a 
significant contribution to the cross section. Thus, the 



TABLE II: Enhancement and suppression of transmission co- 
efficients relative to an equal parity distribution for direct 
capture on a target state with parity 7r^ and ending in a final 
state with parity tv b . The factors JF T and T are defined in 
Eqs. (|32[) and (|33[) . respectively. 



TV 


Allowed Transition 


Factor 


K A = n B 


£a — 1,3,5,... 




Lb 1 V 


i a = 0,2,4,... 





ratio % to be inserted in Eqs. (|32j) and ((33)) becomes 

«/,p-wave 

X = u->i/.s-wave ' (37) 

^DC 

with x < 1 because of the higher s-wave cross section 
(cr s oc 1/y/E) at small energy. Because of the limited 
availability of alternative transitions, contrary to the 
Hauser-Feshbach case the suppression or enhancement 
cannot be compensated in a direct reaction and the fac- 
tors and will have their full effect on each tran- 
sition. For example, \ = 0.18 for capture of a 100 keV 
neutron on the + ground state of 140 Sn and ending in a 
final 3/2 state with an energy release of 1 MeV (see also 
|40|). The inclusion of thermally populated target states 
will slightly reduce the effect by making available states 
with different spin but same parity because of a predom- 
inant parity at low excitation energy. This will have a 
smaller impact than in the Hauser-Feshbach model since 
low partial waves dominate the cross section, there is no 
moderating effect of a total width, and generally fewer 
possible transitions. 

The application rule for J 7 ' and for direct capture 
is shown in Tab. [Ill For the case of the final nucleus in a 
capture reaction having opposite parity than the target 
ground state, the capture rate will be enhanced relative 
to a rate calculated using a level density with equally 
distributed parities because of the predominance of states 
with ground state parity at low excitation energies. This 
may only be altered if an overshooting of the opposite 
parity occurs, as discussed in Sec. Ill Dl However, direct 
capture is important only at such low separation energies 
that the overshooting regime is not reached. 

3. Discussion 

We have calculated parity distributions for all nuclei 
between proton and neutron dripline from Ne to Bi and 
recomputed the Hauser-Feshbach astrophysical reaction 
rates given in [321.I33T] with the same parameters as quoted 
there (set FRDM) . Selected examples were already shown 
in Figs. [T2TFT51 In Sec. IIII C"T1 we discussed considerable 
enhancement or suppression for certain transitions. In 
the astrophysical context, however, the impact of the new 
parity ratio remains limited. This is due to several rea- 
sons. 
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FIG. 16: Stellar reaction rate ratios for 94 Nb(p,7) 95 Mo as 
function of stellar temperature T* (in units of 10 9 K) ; plotted 
is the ratio of the stellar rate ro calculated with the parity 
distribution derived in this work with the standard rate n 
from 32, 33] . Further plotted is the ratio of the same rates but 
without inclusion of experimentally known low-lying levels 
(see text). 
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FIG. 17: Same as Fig.[l6]but for the reaction 95 Zr(p,7) 96 Nb 



Firstly, the impact of a change in the parity distribu- 
tion is higher for reactions with large Q values because 
a larger number of transitions from the compound nu- 
cleus is energetically accessible. For example, neutron 
captures will be altered more strongly on the proton-rich 
side whereas capture on neutron-rich targets - provided 
the compound reaction mechanism is dominating - will 
be less affected. A similar observation can be made for 
proton capture. Thus, the proton captures occuring in 
the rp proce ss lia. [4611 and the neutron captures in the r 
process [13, l48l l49l l50l| show changes of several percent 
relative to the rates computed with an equal distribution 
of parities. 

Moreover, for small capture Q values, the statistical 
model will not be applicable and resonant or direct cap- 
ture will dominate. We expect a larger impact of the 



parity distribution in the averaged direct capture model 
but further conclusions have to await the recalculation 
of the direct rates. Even with changed rates, however, it 
has to be kept in mind that the processes occurring close 
to the driplincs are often in equilibrium, where individual 
rates are not important to determine the produced nu- 
clear abundances [5l| . Rates will only become important 
in the freeze-out phases when the nucleosynthesis path 
moves closer to stability [H, H^, [5l| . In these phases the 
compound mechanism will set in again. 

Close to the valley of stability level schemes are known 
to comparatively high excitation energies and are used 
in the calculations of the relevant nuclear cross sections. 
Usually, up to 20 levels are included whenever known 
[32l [33| . This will further reduce the effect of the applied 
parity distribution because an effective, non-equilibrated 
parity distribution is already in use for the lowest ex- 
citation energies. Figs. fT2T[T5l show a comparison of 
rates with and without included experimental level in- 
formation. It should be noted that the slight change 
in 67 Se(n,7) is not due to neglected levels in 68 Se but 
rather due to the modified transmission coefficients in 
the competing channels, i.e. proton- and a emission. In 
the example shown in Figs. [TBI and [TTl 20 low- lying lev- 
els have been considered in all channels in the standard 
calculation. Removing them, the impact of the parity 
distribution when using pure level densities can be seen. 
Generally, only few low-lying states have to be included 
because the transitions with the highest relative energy, 
leading to the lowest states, dominate in all channels. 
Due to the inclusion of experimental information in the 
standard calculations, the impact of a changed parity dis- 
tribution is also small for the late freeze-out phases and 
the nucleosynthesis phases involving nuclei close to the 
valley of stability. This includes the photodisintegration 
reactions in p-process nucleosynthesis [H, HH l54j . even 
though they occur on a longer timescale than any freeze- 
out process. 

Finally, the most important mechanism suppressing 
the impact of the parity distribution is the necessity to 
include transitions on excited states in the target. In 
an astrophysical plasma, excited states can be thermally 
populated as specified in Eq. ([27)) . Because of the Bohr 
independence hypothesis for compound reactions this al- 
lows for a large number of alternative transitions in case 
one of the parities is suppressed. Small alterations of the 
energy-dependence of the cross sections cannot be seen 
in the astrophysical reaction rates anymore due to the 
folding with the Maxwell-Boltzmann energy distribution 
of the interacting nuclei according to Eq. (|2l)|) . The fact 
that the impact of the parity dependence is strongly sup- 
pressed in stellar rates can be clearly seen in Figs. [T21 fTTl 
The larger the stellar temperature, the smaller the parity 
impact. 

Concluding, inclusion of the parity distribution in this 
work leads to a change in stellar reaction rates of only 
several percent up to about 10% in astrophysically rele- 
vant reactions. A few cases exhibit changes up to a factor 
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of 2 or 1/2, respectively, but are not astrophysically im- 
portant. 

IV. SUMMARY 

We generalized a method to calculate the excitation- 
energy dependent parity distribution in nuclei and used 
it to calculate parity ratios for all nuclei up to Bi. We 
demonstrated the importance of including a sufficiently 
large single-particle model space. In particular, excita- 
tions from core single-particle states can also be impor- 
tant. The parity ratio proved to be non-equilibrated 
up to 5-15 MeV. Interestingly, an overshooting effect 
in the parity ratio was found at major shell closures 



where states with ground-state parity are outnumbered 
by states of opposite parity over a range of several MeV. 

The derived parity distribution was then used to re- 
calculate astrophysical reaction rates. The impact of the 
new description on stellar rates including thermal exci- 
tations of the target proved to be limited in the Hauser- 
Feshbach model although the impact on direct capture 
has to be studied further. 
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